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Abstract 

Astronomical measurements of the Omegas for mass density, cosmological con- 
stant lambda and curvature k are shown to be sufficient to produce a unique and 
detailed cosmological model describing dark energy influences based on the Fried- 
man equations. The equation of state Pressure turns out to be identically zero at 
all epochs as a result of the theory. The partial omega, loa for dark energy, has 
the exact value, minus unity, as a result of the theory and is in exact agreement 
with the astronomer's measured value. Thus this measurement is redundant as 
it does not contribute to the construction of the theory for this model. Rather, 
the value of oja is predicted from the theory. The model has the characteris- 
tic of changing from deceleration to acceleration at exactly half the epoch time 
at which the input measurements are taken. This is a mysterious feature of the 
model for which no explanation has so far been found. An attractive feature of 
the model is that the acceleration change time occurs at a red shift of approxi- 
mately 0.8 as predicted by the dark energy workers. Using a new definition of dark 
energy density it is shown that the contribution of this density to the accelera- 
tion process is via a negative value for the gravitational constant, -G, exactly on a 
par with gravitational mass which occurs via the usual positive value for G. This 
paper also contains an appendix on dark energy dynamics with its own abstract. 
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1 Introduction 

The work to be described in this paper is the generation and study of the theory that 
should go with the experimental workfT] [2] of the astronomers who claim that their 
measurements indicate that the universe expansion is accelerating. Assuming that this 
theory is based on the Friedman equations [TH1 [21]) I have found a solution to those 
equations that seems to be inevitable if the ideas put forward by the dark energy workers 
are to be realized as a consequence of the model. Very briefly their main conclusion from 
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observational astronomy is that a universe expansion process is taking place now that 
changed from deceleration to acceleration at some time t c in the past. This change is 
identified roughly as occurring in association with events at a universe radius r c — r(t c ) 
with an observed red shift here and now in the range, 0.5 < z < 1. I shall use the 
subscript c as indicating the time of change from deceleration to acceleration. They 
assume that the cosmological constant A is positive and deduce from experiment that 
k = 0. 

2 Accelerating Model 

Using Friedman's equations, 

87rGpr 2 /3 = r 2 + (k - Ar 2 /3)c 2 (2.1) 
-8nGPr/c 2 = 2r + f 2 /r + (k/r - Ar)c 2 (2.2) 

in the case k = and a positively valued A we have, 

8nGpr 2 /'i = r 2 - |A|r 2 c 2 /3 (2.3) 
-SnGPr/c 2 = 2r + r 2 jr - \A\rc 2 . (2.4) 

I define a radial length Ra associated with A as 

Ra = 13/AI 1 / 2 (2.5) 

and use Rindler's constant C, 

C* = 8nGpr 3 /3 (2.6) 

to write the first Friedman equation in the form, 

r 2 = C/r + (rc/R A ) 2 . (2.7) 

Having the quantity, Rindler's[16] constant C, in equation (12.71) fixed at a constant value 
turns out to be an important asset when we come to integrate the Friedman equations. 
If as is usually taken to be the case, the mass density function p is taken to be of uniform 
value throughout the universe at any given epoch, t, though having a value that depends 
on epoch through r(t) as p(r(t)) and the quantity 47rr 3 /3 = Vjj(r) is interpreted as the 
total 3-dimcnsional volume of the universe when its radius is r, then the mass of the 
universe will be given by 

Mu(r) - p(r)Vu(r) = 4nr 3 p(r)/3 (2.8) 
and then the constant C takes the form 

C = 2Mu(r)G. (2.9) 
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From equation (|2.9p , it follows while C can be kept constant for variations of the boundary 
radius r of the universe we can incorporate an r dependent gravitation constant, G(r), 
provided an r dependent mass, Mjj(r), is suitably chosen to conform with equation (I2.9[) . 
This point will be returned to when we discuss quantization of the integrated Friedman 
equations. In big bang models in which all the mass of the universe is created at the big 
bang event and thereafter keeps at a constant value Mjj, G must also keep at a constant 
value if the equation (I2.9|) is to be retained for constant C. The equation (|2.9[) can also 
be retained for constant C in cases where the mass of the universe is generated over 
time provided a suitable time dependent G is incorporated to keep C constant over time. 
Most of the work to follow applies to either of these cases so that we do not need to 
specify which form of G is involved. However, we do need to distinguish between the two 
cases if we wish to examine the mass density function p of the universe explicitly. The 
accelerating universe astronomical observational workers [1] give measured values of the 
three fls, and wa to be 

Q Ufi = 87rGpo/(3ff 2 ) = 0.25 (2.10) 
n Ai0 = Ac 2 /(3ff 2 ) = 0.75 (2.11) 
^fc,o = -kc 2 /(rlHl) = 0, k = 0, (2.12) 
= Pa/(c 2 Pa) = -l±«0.3. (2.13) 

Here the value of Hubble's constant will be taken to be 

Ho = 72 Km s^ 1 Mpc" 1 (2.14) 

or in inverse second^] 

H a = 2.333419756287 x 10~ 18 s _1 . (2.15) 

From these values it follows that Rindler's constant C has the form and value 

C = M ,o-ffo r o = 1.361211939757935?^ x 10~ 36 (2.16) 

in terms of the radius now, tq . A and R\ will have the forms and values 

A = A .o3(i/ /c) 2 = 1.363097286965269 x 10~ 52 (2.17) 
R A = |3/A| 1/2 = 1.483532963676604 x 10 26 . (2.18) 

It is convenient to rearrange the first Friedman equation (I2.3[) in the successive forms, 

87rGpr 2 /3 + Ar 2 c 2 /3 = r 2 (2.19) 
87rGpr 2 /3 + H-n:Gp A r 2 /3 = r 2 (2.20) 

1 Many decimal places will be used to keep track of minutely different values that can arise from 
different calculation routes 
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and so identify a mass density p\ which can be used to account for the cosmological 
constant contribution as an additional mass density along with the original p. Thus 



8irGp\ 



Ac 2 



(2.21) 



or 



Pa = Ac 2 /(8ttG). 



(2.22) 



Thus in these terms the first Friedman equation (I2.3P becomes 



8irG(p + p A ) = 3(r/r) 2 = 3H(t) 2 . 



(2.23) 



The second Friedman equation (|2.4p can be written as 



8nGr(~P/c 2 + p A ) = 2f + f 2 jr. 



(2.24) 



The introduction of the additional mass density p\ as in equation (|2.22p to explain 
the mathematical appearance or existence of the cosmological constant A as a physical 
contributor to the theory is the usual approach. However, it is not necessarily the best 
way of physically accounting for the cosmological constant as the resulting equation (12.231) 
has the built in implication of putting p and p\ on a par with respect to the kinematic 
quantity (r/r) 2 = H(t) 2 and this may not be physically very relevant. An alternative 
approach with an alternative density, p A , will be discussed later. 

It follows from the definition (|2 .22[) of p\ that the mass density associated with A 
has the same sign as does A itself which for some theory constructs is taken as positive 
and in other theoretical constructs is taken as negative. The negative mass density case 
does present conceptual difficulties. Using an equation of state involving pressure is also 
conceptually difficult in this context because an equation of state of the form 



for negative A and positive uj\ implies that the pressure Pa is also negative. For positive 
A and negative wa it also implies that the pressure Pa is negative and is the case that 
is used by the astronomers. This negative pressure is usually accommodated with some 
intellectual gymnastics. However, in the case of a partial pressure as in equation (|2.25|) 
it is more easily acceptable^. I show here that such conceptual difficulties are avoided 
in the cosmological model to be examined here. In the next section, I derive a solution 
to the Friedman equations that can accommodate all four of the dark energy researchers 
experimental values. 

2 A confusion as between the total pressure P and the partial pressure in the first version of this 
paper has been rectified in this present version. 



P\ = ^ac 2 pa 



(2.25) 
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3 Dust Model Solution 



From equation (|2.7p . 



(dt/dr) = C,r +{ cr,R A? ^ 

{dt/dr) = ± jcftww^ (3 - 2) 

1 = ± /(C/r + (cr/ J R A )2)i/2- (3 ' 3) 

f ^ r 

± J (C/r) 1 /2(i + ( c / il - A )2 r 3 /c) i/2 ( ■ ) 

± J (C) 1 /2(l + ( c /i ?A )2 r 3 /C) l/2- ( ■ ) 



The change of variable r — » t/ with the standard integration form at equation (|3.9[) gives 
equation (|3.10|l as t in terms of the transformed variable, y. 

r 3 = y 2 (3.6) 

(2/3)d(r 3 / 2 )/dr = r 1 / 2 (3.7) 

r 1/2 dr = (2/3)% (3.8) 

= a - 1 /2l n ( a i/2 2/ + (ay 2 + 1) i/2 )j (39) 



(1 + ay 2 ) 1 / 2 



a = {c/R A ) 2 /C. (3.11) 

Thus the following sequence of steps gives the integral evaluated in terms of the trans- 
formed variable y, an inversion back to the original variable r at (|3.20p and the in- 
troduction of the simplifying function 0± (t) and the constant b to arrive finally at the 
solution 

t = ±{C)- 1/2 {2/3)a- 1/2 ln(a 1 / 2 y+{ay 2 + 1) 1/2 ) (3.12) 
= ±(C)- 1 / 2 (2/3)((c/i? A ) 2 /C)- 1 / 2 
x \n{{{c/R A f/C)^y + {{{c/R A ) 2 /C)y 2 + I) 1 ' 2 ) (3.13) 
±3ct/(2R A ) = Hcy/R A )C-^ 2 + (( C y/i? A ) 2 /C+l) 1/2 )- (3-14) 
{{cy/R A ) 2 /C + l) 1/2 = cxp(±3ct/(2R A ))-(cy/R A )C- 1/2 (3.15) 
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(cy/R A ) 2 /C + 1 = exp(±3ct/(fl A )) -2eM±^t/{2R A ))(cy/R A )C-^ 2 

+ {cylR A ) 2 C-\ (3.16) 
1 = cxp(±3ct/(R A )) 

- 2exp(±3rf/(2i? A ))( C2 //i? A )C" 1/2 . (3.17) 
cy/(i? A C 1/2 ) = ex P( ±3cf /( 2fl A)) ~ exp(T3cf/(2J? A )) 

= sinh(±3ci/(2 J R A )). (3.18) 

cy = R A C 1/2 smh(±3ct/{2R A )) (3.19) 

cr 3/2 = i? A C 1/2 sinh(±3rf/(2i? A )) (3.20) 

r(t) - (i? A /c) 2 / 3 C 1/3 sinh 2/3 (±3ct/(2i? A )). (3.21) 

6 = (i? A /c) 2/3 C 1/3 (3.22) 

0±(t) = ±3ci/(2i? A ) (3.23) 

r(t) = 6sinh 2/3 (6>±(i)). (3.24) 

Equation (|3.2ip or (I3.24[) is a formula for the radius of an expanding accelerating universe 
conforming to all four measurements of the dark energy investigators. 

4 Characteristics of Model 

The velocity, v(t) = r(t) = dr(t)/dt, of expansion as a function of t is given by 

v(t) = (i? A /c) 2 / 3 C 1/3 (2/3)sinh- 1/3 (0±(i)) 

x cosh(0 ± (i))(±3c/(2# A )) (4.1) 

= ±(Cc/i? A ) 1/3 sinh~ 1/3 (0±(i))cosh((9 ± (t)) (4.2) 

= ±(6c/i? A )sinh" 1/3 (6» ± (i))cosh(6' ± (<)). (4.3) 

The acceleration a(t) — v(t) — r(t) is found to be 

a(t) = ((c/i? A ) 4/3 C 1/3 /2)(3sinh 2/3 (0±) 

- cosh 2 (9±{t)) sinir 4/3 (0±(*))) (4.4) 

= 6(c/(i? A )) 2 sinh 2/3 (6» ± (t))(3-coth 2 (6' ± (t)))/2. (4.5) 

From equation (|3.21[) and equation (|4.2p . Hubble's constant is given as a function of t as 

= ±(Cc/R A )^sinh- 1 / 3 (6 ± (t))co S H9 ± (t)) 

(i? A /c) 2 / 3 C 1 /3 S inh 2/3 (0±(t)) 
= {c/R A )coth(±3ct/{2R A )). (4.7) 
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The measured value of H{t) is given by equation (|2.15l) . Thus a consistent value for time 
now, to, is the time solution of (14.81) at (|4.9p with a numerical value in seconds at (|4. 10|) 

(c/R A ) coth(±3ci/(2i? A )) = 2.333419756287 x 10~ 18 (4.8) 
t = (±2R A )/3c) cotlT 1 ((R A /c) x 2.333419756287 x 10~ 18 ) (4.9) 
= 4.34467334479058 x 10 17 s, (4.10) 

using the value of R\ from equation (|2.18|) . From the formula for the acceleration (I4.3[) . 
it can be seen that there is a time t c at which the acceleration changes from negative to 
positive given by the time solution of (|4.11l) at (|4.12p 

cosh 2 (6>(i))sinh- 4/3 ((9(i)) = 3 sinh 2/3 (0(t)) (4.11) 
t c = (2i? A /(3c))coth _1 (3 1/2 ) = 2.172336672394881 x 10 17 . (4.12) 

From equations (|4.17[) and p.21[) we can find the value of r = r(t ), the radius of the 
universe formally at time-now as 

r = (i? A /c) 2/3 C 1/3 sinh 2 / 3 (±3ci /(2i? A )). (4.13) 

The evaluation of which requires that we have a value for Rindler's constant C . However, 
an inspection of the definition for this constant equation (|2.6p shows that to know this 
constant we require knowing Hq for which I have assumed an input value at equation 
(|2.15[) but it also requires knowing the present radius of the universe ro for which we have 
not got a measured value. There is an alternative approach to finding the value for the 
present time, to, and radius, ro = r(to), of the universe and thus getting the appropriate 
value for Rindler's constant. From equations (|2.6[) and (|2.10p . we note that Rindler's 
constant can be expressed in the form 

C* = ^Mi^V 3 , (4.14) 

where I am now using the dagger superscript to denote the present day radius of the 
universe and the present day value for Hubble's constant. The dagger notation is to 
emphasise the alternative route used to calculate quantities such as r T . From equation 
(|3~21"j) for r(t) it follows that 

r t3 = (i? A /c) 2 C , sinh 2 (±3ci t /(2i? A )). (4.15) 

The value for Rindler's constant, C, from equation (|4. 14[) can be substituted into equation 
(|3T5)) to give 

1 = 0.25(i? A /c) 2 ff t2 sinh 2 (3ctV(2i? A )) (4.16) 
having cancelled the r^ 3 = r 3 (t') that occurs on both sides of the equation. Thus 

± 1 = Q.h{R K /c)H^ sinh(3rfV(2i? A )). (4.17) 
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This can be solved for tf. Taking the plus sign with unity and using equation (|4.7[) at 
time now as in (14.181) . 

= H(P) = (c/R A ) coth(±3c*V(2flA)). (4.18) 

It follows that the value of ft is given by the time solution of (14. 19|) as t< at (|4.20l) . 

2 = cosh(6»(t t )) (4.19) 
t f = (2R A /3c) cosh" 1 (2) = 4.344673344789258 x 10 17 s. (4.20) 

Remarkably close to the value obtained for the original to displayed again at equation 

<o = 4.34467334479058 x 10 17 s. (4.21) 

However, the value of v was obtained without directly using the numerical value of 
H at equation (|4.6|l although it does indirectly depend on H through the value of R\ 
so that it can be used to check the theoretical numerical value of Hubble's constant 
W = H(t>) using the theoretical formula for H(t). Thus the time v is greater than t c , 
the acceleration change time, by an amount 

t + - t c = 2.172336672394881 x 10 17 s = 6.888434400035 x 10 9 y. (4.22) 

The v notation from now on will be used for time-now rather than the original to but 
this implies that the value of Hubble's constant should be given by formula (|4.7|) as 

ff(ft) = (c/R A ) coth(±3c^/(2i?A)) = 2.333419756287 x 10~ 18 s" 1 . (4.23) 

Comparing this with the experimental value 2.333419756287 x 10 -18 s _1 at equation 
(|2.15[) . we see that the value coming from this theory via the indirect route through the 
value of R A is exactly the experimental value given at (|2.15|) . However, although we have 
found v at formula (|4.22p . we still have not found the value of r' = r(v) that is necessary 
to find Rindler's constant C and the constant b. It seems that mathematically C is to 
be taken as arbitrary or alternatively is to be taken as arbitrary or either C or is to 
be obtained from experiment or further some other theoretical consideration needs to be 
used to obtain one or other of these two key constants. There are values ascribed to 
from astronomical observation based on a various speculative extrapolations and which 
are therefore not greatly reliable see Gravitation[T7] page 738, box(27.4). I have here 
decided to use a value for r* that comes from an assumption about the dependence of 
the gravitation constant on v or on ct^ . This assumption comes from a suggestion about 
a formula for G in terms of other physical constants of great numerical accuracy noticed 
by Ross McPherson, see references [HI [14]. His original suggestion had dimensions 
different from G. A generalised version of McPherson's suggestion but which has the 
usual dimensions associated with the gravitation constant G is as follows 

G = h 2 /(m 2 m e ct ), (4.24) 



8 



where m p and m e are the rest masses of the proton and electron respectively. A formula 
similar to (|4.24[) for the gravitation constant was suggested by P. A. Dirac[13] in the 
thirties. There are a number of ways in which this formula can be interpreted. Taking to 
to have one of values associated with the present age of the universe gives a very accurate 
value for G[TT1 [TBI HO]- Otherwise ct can be given a value associated with the present 
radius of the universe r* multiplied by cos(xg), the numerical constant, cos(xg) ~ 1 
coming from gravitational coupling, which will also give a value for G. I have used this 
last option to develop a quantum theory for gravity[23 1121 IS] that successfully gives 
accurate formulae conformable with Dirac's large number hypothesis. This gives me 
some confidence in now inverting the formula (j4.24[) and using the measured value for G 
to supply the missing numerical value for as 

rt = h 2 /{m 2 p m e G) = 6.539532681821 x 10 25 , (4.25) 

where the numerical constant, cos(xg) ~ 1, has been omitted as its values is very close to 
unity at nearly all times. The situation when it could make some significant difference, 
at the start of expansion of the universe when t w +0, will be discussed in a future paper. 
The issues of reformulating quantum theory so that it is consistent with the theory of 
stochastic processes, the quantization of gravity theory and the consequent production 
of a quantized cosmology is discussed in references [H [6j [7J HI [10]. The use of 
the formula (|4.24[) for G in terms of quantum quantities in this work represents a weak 
partial quantization of the structure. A more complete quantization will be discussed 
in a later paper. This can now be used to complete the numerical computation for the 
present theory by using the value for in equation (|4.25p to give us Rindler's constant 
C through equation (|2.16|) which is repeated below with r' replacing tq at (|4.26|) and its 
numerical value is given at (|4.27|) . The value of the constant b is then deduced at (|4.28|) . 

C = ft M ,oH^ 3 = 1.361211939757935r t3 x 1(T 36 (4.26) 
= 1.361211939757935r t3 x 1(T 36 = 3.806871984611 x 10 41 (4.27) 
b = (i? A /c) 2/3 C 1/3 = 4.53 4 2 5 8 7 1 39 2 5 x 10 25 . (4.28) 

I finish this discussion with a few remarks about the parametric values that arise in this 
theory. Firstly we can calculate the red shift, z, that would apply to light emitted at the 
acceleration change time t c from the universe boundary r(t c ) to arrive at a present day 
observer as 

z = r(J)/r{t c ) - 1 « 0.81712. (4.29) 

This certainly strengthens the theory-observation cohesion in that it is a theory result 
within the limits suggested by the dark energy workers. An interesting and curious direct 
result from this structure is that the ratio of time now to the time of acceleration change 
as calculated from their respective values is exactly 2, 

J/t c = 2. (4.30) 
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One might have noticed the value of this ratio from equations (|4. 12[) and (|4.20|) , if one 
had at some time earlier come across the exact numerical relation (|4.36l) . 



coth 2 (6>(t c )) = 3, (4.31) 

cosh(6»(t t )) = 2, (4.32) 

sinh(6»(t t )) = 3 1/2 , (4.33) 

H(t c ) = 3 1/2 c/R A (4.34) 

= ccoth(6»(t t ))/i?A (4.35) 

cosh _1 (2) = 2coth~ 1 (3 1/2 ). (4.36) 

H(t) = H(t c )coth(9(t))/3 1/2 (4.37) 

H(tf) = 2H{t c )/3. (4.38) 



The first equation above (|4.3ip essentially defines the time t c of change of acceleration 
from negative to positive and comes from equation (|4.12[) . The second equation above 
comes from equation f!4. 17|) or (|4.19l) and defines the time-now , t\ when the input mea- 
surements were made. The fourth equation (|4.34[) comes from equations f|4. 18[) and 
(|4.31|) . I have not been able to decide whether or not the result equation (|4.30|) is some 
remarkable coincidence in values of the present day v and the past t c or some deep 
indication of a need to have more than just the past epoch value t c to be able to talk 
about a definite location in the time history of an event such as a universally distributed 
change from deceleration to acceleration. However, equation (|4.37l) does clearly show 
up the background mathematical theory structure of the relation between t c and t at 
equation (|4.30p in relation to the time evolution of the system and gives equation (|4.38l) 
the Hubble equivalent of (I4.30[) when t = v but, I have to admit, that I have not yet 
deciphered the physical significance of this mathematics as expressed by (|4.30|) or (I4.38|) . 
It would not seem so remarkable if they, (I4.30P and (|4.38[) . were just approximate results 
but in the construction of this model they are numerically exact and unexpected. 

The formula for the pressure from the Friedman equation (|2.4[) is, using the definitions 



for H and R A at equations (|4T6]) and (|2"TTg|i in (j4~4"T]) . 

P = -(c 2 /(8nGr))(2r + r 2 /r-\A\rc 2 ) (4.39) 

= -(c 2 /(87rG))(2r/r + r 2 /r 2 -\A\c 2 ) (4.40) 

= (-c 2 /(8nG))(2r/r + H(t) 2 -3(c/R A ) 2 ). (4.41) 



From equations ((3~2"4"1) . (|43|) and (|4~5j) we have 

2r(t)/r(t) = (c/ R A ) 2 (3 -coth 2 (9+(t))) (4.42) 

H(t) 2 = {c/R A ) 2 coth 2 (9 +(t)) (4.43) 

2r(t)/r(t)+H(t) 2 = 3(c/R A ) 2 . (4.44) 



Substuting these formulae into the form for pressure equation (I4.41[) . we find complete 
cancellation on the right hand side to give an identically equivalent to zero total pressure 
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at all times, 



P(t) = (3c 2 /(8ttG))(~(c/R a ) 2 + (c/R a ) 2 ) (4.45) 
P{t) = 0^w = 0. (4.46) 

Thus this system rigorously describes a dust universe. From equation (14.41| or (|4.45l) we 
see that the total pressure P can be expressed as the sum of the two partial pressures 
Pa and Pm, the first of negative sign and the second of positive sign, as 

Pa = - (3c 2 /(8^G))(c/i? A ) 2 < (4.47) 
Pm = (3c 2 /(8irG))(c/R A ) 2 > (4.48) 
P = Pm+Pk- (4.49) 

The choice of which pressure is associated with which part of the total pressure being 
determined by the sign of that part so that normal gravitational mass gives negative 
acceleration or positive pressure, that is gravitational attraction. The lo a that goes with 
the A equation of state involves negative pressure or positive acceleration is thus given 
by 

- Pa/(c 2 Pa) = -(3c 2 /(&7rG))(c/R A ) 2 /(c 2 p A ) (4.50) 
= -1 (4.51) 

an exact value from theory for uj a at the centre of the measurement range at equation 
(|2 . 1 3|) . A simple interpretation of the force structure that is the cause of the acceleration 
(|4.5|) that is operative in this model is obtained if we introduce a mass density to 
account for the cosmological constant A as 

p\ = Ac 2 /(4vrG) = 2 PA . (4.52) 

That is to say its values is twice the value of the usual density function at (|2.22p . The 
Friedman equation (|2.4p repeated below at equation (|4.53[) with P = 0, can then be 
expressed as equation (|4.54[) through to equation (I4.59[) . 






= 2f + r 2 /r - Arc 2 . 


(4.53) 


r 


= Arc 2 /2-r 2 /(2r). 


(4.54) 




= Krc 2 /2 - {AirGpr + Ac 2 r/2)/3 


(4.55) 




= (Ac 2 - 47rGp)r/3 


(4.56) 




= 4^rG(pi-p)/3 


(4.57) 




= 4^r 3 G(pi - p)/(3r 2 ) 


(4.58) 




= M\Glr 2 - MjjG/r 2 


(4.59) 


M\ 


- ^r 3 p{/3 


(4.60) 


M v 


= 47rr 3 p/3 


(4.61) 



11 



Thus associating the density function p' A and consequent total mass M\ within universe 
with the cosmological constant A we get the very transparent formula for the dynamics 
under gravity, the acceleration that any particle would experience at the boundary of the 
universe being given by equation (|4.59p . It tells us clearly that the normal mass density 
p gravitationally causes the usual decelerating attraction of this particle to within the 
universe body whilst the dark energy mass density causes an acceleration that repels this 
particle towards the outside of the universe. The theoretical structure described here 
involves this as a naturally occurring effect which is built into this model. Consequently 
any need to introduce negative mass densities to describe the dark energy contribution 
together with the conceptually difficult concept of negative total pressure are removed 
from contention. 



5 Conclusions 

It has been shown that a subset of the measurements by the dark energy workers given 
at equations (|2.10[) . (|2.11l) and (|2.12l) together with an input value for Hubble's constant 
now and the assumption that Rindler's constant, C, is an absolute constant lead to 
a unique solution of the Friedman equations. The fourth measurement from the dark 
energy observations at equation (|2.13[) is not necessary for finding this solution. In fact, 
we have found that this last measurement with the partial wa value, lo\ = — 1, is an 
exact result derivable from the unique model from the first three measurements. Thus 
the model derived from general relativity via Friedman's equations is conformable to all 
four of the dark energy measurements (|2.10l) . (|2.11[) . (|2. 12|) and (|2.13j) . 
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Appendix Abstract 

It is shown that negatively gravitating particles can consistently be con- 
sidered to exist and interact with normal positively gravitating particles in 
the contexts of general relativity and classical Newtonian gravitational the- 
ory. This issue arises from the discovery of dark energy which is considered 



12 



to be causing an acceleration of the expansion of the universe. The issue is, 
can this dark energy occur in particulate form? A related issue was studied 
in the fifties by Herman Bondi, (|41j) when he investigated the possible exis- 
tence of negative mass particles in the general relativity context, long before 
dark energy appeared on the scene. He came to a paradoxical conclusion 
that seemed to rule out the actual physical existence of negatively gravitat- 
ing particles. This paradox does not occur in this work because only positive 
mass particle are involved whatever their gravitational character may be. 
The structure of the differential equations that would apply in the case of a 
binary pair of opposite gravitational character components are used to show 
and explain how they can become consistent in general relativity or classical 
gravitation theory. This involves explaining a non-obvious relation between 
the principle of equivalence and Newton's Action Equals Reaction principle. 
A path structure for a mixed mass binary pair is set up which satisfies the 
equations of motion and does not have paradoxical properties. The force 
structure of the system is checked with a known classical dynamical test 
for the force per unit mass involved in the component particles motions. 
This test is used to demonstrates that the basic assumptions of this theory 
are incorporated into the consequential orbital structure. An alternative to 
the Action Equals Reaction principle more appropriate to the astronomical 
situation is suggested. An animation using Mathematica has been derived, 
and is available, and shows how a mixed gravity binary pair move under 
their mutual gravitational action. 

6 Appendix Introduction 

The work to be described in this paper is an application of the cosmo- 
logical model introduced in the papers A Dust Universe Solution to the 
Dark Energy Problem [23], Existence of Negative Gravity Material. Iden- 
tification of Dark Energy [21] and Thermodynamics of a Dust Universe 
[32] . The negatively gravitationally characterised mass density involved 
from which such hypothetical particles might be formed has the value, p^, 
which is twice the Einstein dark energy density, Ac 2 /(8ttG). The issue 
to be examined in this paper is the form that the dark energy negatively 
gravitating material is likely to take. I shall assume that it is in the form 
of material particles like the normal positively gravitating particles and 
consider what consequences this has for the dynamics of this exotic ma- 
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terial in relation to the normal gravitating material. The previous work 
on this topic does lead towards making this particulate assumption about 
the form of dark energy by showing that with the new definition, p' A , for 
dark energy density it comes physically onto a par with positively grav- 
itating material simply by its gravitating strength being denoted by -G 
times mass in place of +G time mass as is usually the case. All of this 
work and its applications has its origin in the studies of Einstein's gen- 
eral relativity in the Friedman equations context to be found in references 
W, |22], [21], [20], [19], [E], [1], [23]) and similarly motivated work in refer- 
ences ([10], [9], [8], [7], [5]) and ( [12] , [13] , [H] , [IS] , [7] , [25] , [3] ) . The applications 
can be found in ( [23J , [24] , [32] , [36] , [31] [10] ) . Other useful sources of informa- 
tion are ([UJ,[3],[30j,[27j,[29],[28]) with the measurement essentials coming 
from references ([I],[2],[TT],[37]). Further references will be mentioned as 
necessary. 



7 Positive, Negative Gravitating Particle Dynamics 

The term pure particle will be used with an extended sense of meaning either 
astronomical accumulations of mass with order of the amount of mass that 
might be found within a planet, a star or a galaxy or quantum elementary 
particle and accumulated structures formed from elementary particles such 
as atoms or molecules. In other words, a pure particle can be mass-wise any 
physical object formed from a mass accumulation with some restrictions. I 
shall assume that particles can have the two possible states of either being 
gravitationally positive, m + , and consequently causing acceleration towards 
itself on all other particles in its vicinity or being gravitationally negative, 
m~ , and causing an acceleration away from itself on all other particles in its 
vicinity. All the particles involved are assumed to have positive mass and 
the restriction pure means that all such particles, whatever mass size will 
each have a definite gravitational coupling value, ±G, G being Newton's 
gravitational constant. That is to say, for a massive particle of mass rri^ 
gravitating strength equals ±Gm . This is a provisional restriction that 
will certainly need to be lifted that I introduce to avoid mixed character 
particles that would involve changed values for G. Further restrictions on 
this very wide definition for possible masses will be discussed later. I shall 
designate the gravitationally characterised positive particles as being of red 
colour and gravitationally characterised negative particles as being of green 
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colour for ease of reference and diagrammatic use and further represent the 
red particles with a + superscript and the green particles with a — su- 
perscript. I emphasise that both types of particle have positive rest mass. 
Thus a red particle induces, at distant points, accelerations towards itself 
in all other particles in its vicinity. A green particle induces, at distant 
points, accelerations away from itself in all other particles in its vicinity. 
This assumption about gravitational characterisation is very different from 
the electric characterisation for charged particles with charges ±e where 
oppositely characterised particles mutually attract and same characterised 
charges mutually repel, the distinction emphasised in the all for the gravi- 
tational case. This may seem to the reader to invite a paradox but it will 
be explained why this paradox does not arise. The classical Newtonian 
equations that describes how a red particle behaves dynamically under the 
gravity field of another red particle are 

f = m\d 2 vjdt 2 = -Gm+m+r/r 2 (7.1) 

d 2 r/dt 2 = (r - r6 2 )r + (rO + 2r0)t (7.2) 

a rtR = r-d 2 r/dt 2 = (f - r6 2 ) = -Gm+ / V 2 : . (7.3) 

a e , R = i-d 2 r/dt 2 = (r6 + 2f6) = r- 1 (d(r 2 6)/dt) = 0. (7.4) 

The force acting is the gravitational force originating from the positively 
gravitating particle mj, regarded as fixed at its specific origin of coordi- 
nates, is directed towards it. That is to say, m\ is at the tail of the vector 
r and m\ is at the sharp end of r and usually in motion. These equations 
contain much information about classical Newtonian gravity. Very striking 
is the fact that the two components of acceleration a r ^ and ag t R of the 
particle, mf, do not involve the mass of this particle, in fact it could be red 
or green and in either case the equation would be unchanged. The mass of 
the red particle is involved in the r component through the influence of 
the gravitational field. Thus the accelerating effect on any particle of any 
mass or any other characteristic under the influence of a red gravitating 
source is the same for all particles passing through the same position dis- 
tant from m\ at any time. Thus from classical dynamical theory, we obtain 
the essence of Einstein's 'principle of equivalence. The main conclusion, I 
wish to emphasise is that the acceleration, a r> R, is always negative with di- 
rection, — r that is towards rn^, for any particle experiencing the effect of a 
positively characterised particle such as because G and m^/r 2 are both 
always positive. That is the induced acceleration is towards a positively 
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gravitating source, m\ whatever the gravitational character the subjected 
particle, mf , may have, 

ar,R= ¥ -^ = (r-r6 2 )<0. (7.5) 

Let us now consider the case of a green particle fixed at the origin of co- 
ordinates being the source of the gravitational force on a red particle. For 
this case equation 17.31 has to be replaced with the negative — G becoming 
+G as in the next equation. 

f - r 

a rG = — = (f-rO 2 ) = +Gm^/r 2 . (7.6) 
m{ 

This then implies also from equation (17. 3p 

a r ,G= t -^ = (r-r6 2 )>0, (7.7) 
rn{ 

as all rest masses and G/r 2 are always positive. Thus a negatively gravi- 
tating mass source particle induces acceleration away from itself for all par- 
ticles in its vicinity regardless of any other characteristics they may have. 
Particle, , in equation 17. 61 differs from the positively gravitationally char- 
acterised particle m\ in equation 17.31 in that it is negatively gravitationally 
characterised. It follows that for a binary system composed of a positive 
gravitating mass and a negative gravitating mass we need two equations of 
motion such as, 

f x = m fd 2 r 1 /dt 2 = +Gm J [m 2 Y l /r 2 l (7.8) 
d 2 ri /dt 2 = (h - rjDh + (n0*i + 2f 1 1 )t 1 (7.9) 

a n , G = t ^ = (r 1 -r 1 e 2 1 ) = +Gm 2 /r 2 1 >0. (7.10) 



a 9l , G = fj -^ = (r 1 9 1 + 2f 1 6i) = Q, (7.11) 
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f 2 = m 2 d 2 r 2 /dt 2 = — Gm\m 2 r 2 /r 2 . (7-12) 
d 2 r 2 /dt 2 = (r 2 - r 2 6 2 2 )r 2 + (r 2 2 + 2r 2 2 )i 2 (7.13) 

a r2 ,R = ^ = (f 2 -4 2 ) = -Gm+/r 2 2 <0. (7.14) 
m 2 

ae 2 , R = = (r 2 6 2 + 2r 2 6 2 ) = 0, (7.15) 

The two sets of equations above refer to a system of two particles one, 
m{, with a positive gravitating characteristic and one, m 2 , with a negative 
gravitating characteristic, both have positive rest mass. In the first set 
is fixed at the origin of coordinates and is the source of the gravitational 
field that determines the dynamics of mf which is found at vector position 
ri at time t relative to the position of m 2 . In the second set m\ is fixed 
at the origin of coordinates and is the source of the gravitational field that 
determines the dynamics of m 2 which is found at vector position r 2 at time t 
relative to the position of mf. Thus there are two space frames of reference 
involved with their origins separated by the relative position vector, r = 
i"i — r 2 , of the particles. The key quantities that determine the kinematic and 
dynamic behaviour of the total system are the four component accelerations 
«r 2 ,-R> a n,R, a r 2 ,G, &ri,G but being measured in different frames of reference 
it is difficult to see how their effects are to be combined. This is necessary 
if we are to understand how the particles interact. However, Inspection of 
the two basic equations of motion (17 ,8p and (I7.12p . repeated below, which 
somehow or other are to hold as a simultaneous pair for the two particle 



system under examination, 

fj = m+dhjdt 2 = +Gmfm 2 r 1 /r 2 1 (7.16) 

f 2 = m 2 d 2 r 2 /dt 2 = —Gmfm 2 r 2 /r 2 1 (7-17) 

h =- r 2 (7.18) 



the first impression is that they are incompatible as they seem to defy the 
Newtonian law that action and reaction should be equal. The first equation, 
gives the force from m 2 acting on m\ . The second equation, gives the force 
from mi acting on m 2 . Usually, if the action is taken to be f~i, the reaction 
is regarded as being f 2 , then, f x = — f 2 , but here because of (" 17. 1 8[) it follows 
that fi = f 2 . This is a well know paradox arising when trying to work 
with negatively gravitating particles. Subtle, you may say but nevertheless 
a source of great consternation in this area of work. It is clear from these 
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equations that two negatively gravitating particles in interaction would also 
give the result fi = — f 2 , implying that their interaction satisfies Newton's 
third law so that we only have to consider the two situations of positively 
gravitating pairs and oppositely gravitating pairs to analyse and explain 
this situation. Thus we now have to consider the way this difficulty can be 
circumvented. We can get some direction in this investigation by considering 
how the similar well known problem in the study of the orbiting positive 
gravitating particles of a binary system in astronomy is handled. There are 
also two equations of motion as in (I7.16[ I7.17p . 

f x = mfdhjdt 2 = -Gmfm+ri/rl (7.19) 
f 2 = m +d 2 Y 2 /dt 2 = -Gmfmth/rZ. (7.20) 

Clearly there are no negatively gravitating particles involved here and here 
we do have action and reaction equal because f 1 = — r 2 which implies fi = 
— f 2 . The usual procedure is to get the two equations combined by defining 
the centre of mass, r cm , frame in place of the two obviously different frames 
used above as follows: 

mtvi + mfr 2 

T = 

m{ + 

In the absence of any external forces acting on a binary system, it is as- 
sumed that the center of mass can acquire no acceleration. Thus we can, 
by differentiating twice with respect to t the centre of mass vector, deduce 
that 

= d 2 r cm /dt 2 = m^d 2 ri/dt 2 + m 2 d 2 Y 2 /dt 2 , 

as the centre of mass denominator factors out. It follows that 

= d 2 r cm /dt 2 = mfdhjdt 2 + m+d 2 r 2 /dt 2 = fi + f 2 = 0. (7.23) 

Thus the relations above become identities because the usual action equals 
reaction condition holds for the classical binary particle pair case. Clearly 
also, this condition as it stands does not hold for the mixed character particle 
pair case and it is fairly obvious why this is so. The centre of gravity has to 
replace the centre of mass in the mixed mass case. However, notably, in the 
classical situation the centre of mass and the centre of gravity of the systems 
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(7.21) 



(7.22) 



coincide. The centre of gravity in the mixed mass case can be defined as 

G+m^Yi + G-m 2 r 2 _ mj^ - m 2 "r 2 
Lx + m{ + G_m 2 m{ — m 2 

G + = +G (7.25) 

G_ = -G. (7.26) 

In the mixed mass case, differentiating r cg twice with respect to t and as- 
suming that the centre of gravity cannot acquire any acceleration in the 
absence of external forces gives 

= d 2 r cg /dt 2 = m+d^/dt 2 - m^d 2 r 2 /dt 2 = f x - f 2 = 0, (7.27) 

because in this case, the condition, fx — f 2 — 0, contradicting action and 
reaction are equal holds. In the light of the form these structures take, 
it is perhaps tempting to conclude that all that need be done to get a 
consistent theory is to redefine the equality of action and reaction as the 
magnitude of action and reaction forces are always equal. I am reluctant 
to do this for reasons to be discussed later. Using the basic two types of 
structure either centre of mass or centre of gravity orientated together with 
the separation vector of the two particles concerned r(t) = ri(t) — r 2 (t) the 
position vectors for the two cases (I7.23P and (17.271) can respectively be each 
expressed in terms of r(t) in two equations, 



r 2 



ri - r cg - 



r2 - r cg - 



mj + 


(7.28) 


mj + 


(7.29) 


+ rrt) 
mj — m 2 


(7.30) 


mt . . 

+ -r * • 
mj — m 2 


(7.31) 



If we now differentiate these four equations through twice with respect to 
time assuming that the acceleration of the centre of mass or the centre of 
gravity is zero and then multiply through by the appropriate gravitationally 
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subjected particle we get 

mfd^/dt 2 = + ^ m ^ rf 2 r(0M 2 = fi = -f2 (7.32) 

m 2 + rf 2 r 2 /rft 2 = - m / m \ d 2 v{t)/dt 2 = h (7.33) 

m7/ + mjj" 

mtd^/dt 2 = -^$^-d 2 r(t)/dt 2 = f 1 = +f 2 (7.34) 
m7/ — m 2 

m 2 d 2 v 2 /dt 2 = - m2m+1 d 2 r(t)/dt 2 = f 2 . (7.35) 
m~[ — m 2 

The last four equations can be used to summarise results so far. The first 
two equations above are well known results from Newtonian gravitational 
theory and are used to find the astronomical orbits of binary star systems. 
The last two of the four are new results that describe the dynamics of bi- 
nary mass systems composed of one positively gravitating mass together 
with one negatively gravitating mass and it has been shown that the new 
system seems to defy Newton's action and reaction law. This being the puz- 
zling situation that the normally gravitating mass attracts the negatively 
gravitating mass whilst the negatively gravitating mass repels the normally 
gravitating mass. It has long been thought that this is a non reconcilable 
paradox that excludes the existence of negatively gravitating particles from 
the physical arena. An important point I wish to emphasise is that the 
negatively gravitating mass in this theory does not have negative mass. The 
gravitational negativity is an intrinsic property of its positive mass struc- 
ture just as the negative charge of an electron is an intrinsic property of 
its positive mass structure. From the equations above the two component 
systems reduce to a single component equation in r(t) for the normal com- 
ponent pair and the mixed component pair. We note also that the usual 
reduced mass that arises in the binary system, miG + m 2 /(G + m 1 + G + m 2 ), 
in the mixed mass system is replace by miG m 2 / \G +mi + G_m 2 ) and this 
will change sign if the mass character of the masses are interchanged. The 
form of the reduced mass for the mixed mass system also exposes a difficulty 
with the single equation of motion that arises for the mixed mass system. 
If the two mass components and m 2 in such a system are equal their 
reduced mass becomes infinite. Thus rendering the reduced mass equation 
of state unusable. Consequently a binary mixed gravitating pair with equal 
masses must be excluded from discussion. I am sure that this difficulty has 
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some fundamental significance but I do not know what it is. In the next 
section, I find a full solution for the paths of a mixed mass pair under their 
mutual gravitational interaction. 



8 Mixed Mass Pair Paths 



In this section, I shall derive a complete solution to the dynamical problem 
of a positively gravitationally characterised mass and a negatively charac- 
terised mass moving in conjunction under their mutual gravitational inter- 
action. In the previous section, it was shown that such a system seems 
to violate the usual version of Newton's law of action and reaction being 
equal. I shall work with a special case and show the solution does not in- 
volve bizarre features such as one particle chasing another to infinity which 
arose in the Bondi analysis. The production of this solution, I regard as 
something like a mathematics existence theorem, here showing that nega- 
tively gravitation particles can occur in nature in interaction with positively 
gravitating particles with unambiguous orbits, with a rational explanation 
and not contradicting general relativity or Newtonian dynamic. Repulsive 
inverse square law force is well known in the electromagnetic theory in the 
context of the electron and positron interaction for example. The hyper- 
bola is particularly interesting in the inverse square law context as it has 
two branches each with it own focus and with reference to one or other of 
its two foci, the two branches can represent an attracted particle path and a 
repelled particle path. In gravitation, theory the repelled particle path has 
hitherto been regarded as of no interest because gravity has always been 
thought to be only attractive. To give these remarks some mathematical 
basis let us first consider the pedal, (r,p), equations of the two hyperbolic 
branches for a particle moving under inverse square law gravity, 



The first equation above represents the orbit or branch of an hyperbola 
occupied by a particle being attracted to the focus within the orbit. This 
orbit is concave to the active focus. The second equation above represents 
the orbit or branch of the hyperbola occupied by a particle being repelled 
from a focus outside the orbit. This orbit is convex to the active focus. 



{b/pf 



2a/r + 1 
1 - 2a/ r' 
{x/af - {y/bf. 




1 



(8.3) 
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The second equation, (18.21) would usually be rejected in the gravitational 
context. The third equation involves both the branches and is represented 
in a frame of reference with the centre of the hyperbola at the centre of 
coordinates. Even though the receptor particle is negatively gravita- 
tionally characterised, the numerical value of the attractive gravitational 
force from the particle mf, fixed at the local focus, on a particle on the 
first branch above is given by, 

, momtG rruh 2 . . . 

f* = r 2 = -^ d Pl dr i ( 8 - 4 ) 

where h = r 2 8 is the constant arrived at by integrating equation (17.111) or 
(I7.15p . Integrating equating (I8.4p with respect to r, we get 



If we multiply equation (18. 5 j) through with (b/h) 2 , we get 

and comparing (18.61) with (I8.ip . using the relation, b 2 = a 2 (e 2 — 1), between 
b, a and e, the eccentricity, we find, 

° ~ [bj ~ a 2 (e 2 - 1) " h 2 - ffllG/ ° (8 - 7) 

b\ 2 + „ a 2 (e 2 -l)mtG 
' in |' ( . - 



hj 1 h 2 

(8.9) 



(e 2 - \)mjG 

Let us now consider the second branch of the hyperbola given by (I8.2p . 
Firstly, suppose that we are not sure about the force that would have to be 
at the focus of the first branch to control the motion followed by a particle 
on the second path. Thus let us call this unknown or uncertain force 
force G?. We now have to integrate the formula, 

mim-G? rath 2 , . , . . 

h = -^t- = -jTdp/*r> (8-10) 
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where h = r 2 8 is the constant arrived at by integrating equation (17. lip or 
(I7.15p . Integrating the equation (I8.10p with respect to r, we get 

p j r 

If we multiply equation (18.1 ip through with (b/h) 2 , we get 
b\ 2 fb\ 2 2m~G, „fb^ 2 



pj \h J r ~^~^\h) ) 

and comparing (18 . 1 2[) with (18.21) . using the relation, b 2 = a 2 (e 2 — 1), between 
b, a and e, the eccentricity, we find, 

C = (tf-V ^-^' ^c,,* (8,3) 



bj a 2 (e 2 -l) ft 2 
. = -(-j ,n 2 G, = - J (8.14) 

" = - (e^-L-G/ (8 ' 15) 

The two branches refer to the same hyperbola, (18.131) . so that the value for 
a obtained by the two different routes followed above should yield the same 
result. Thus from ( 18. 9 p and ( 18. 15ft it follows that 

& first route -. ^2 ? [q ~[ r\ 
= 1 = T7T- (8.1b) 

" second route Tfl-y Lx + 

Here we seem to have a big problem, because we have been assuming that the 
force induced by a negative gravitating particle at a distant point involves 
G- = — G rather than G + = G and here we see that substituting G_ for 
G? gives the result 

O first route ^ ^2 

(^second route Tfl-y 

This seems to be disastrous because it implies that the two masses have 
to be equal, a situation excluded earlier. However, it can also be taken to 
imply that we require two different sized hyperbolae each with its own a 
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value to get a consistent binary system. Thus we conclude that G? can be 
taken to equal to G_ = — G provided we work with one hyperbola with a 
major axis clr and a second hyperbola with a major axis ac and such that 



^ = 1^ 1. (8.18) 

However, the conclusion that two hyperbolae are necessary for the construc- 
tion of a binary gravitational system is otherwise obvious. It would seem 
impossible to set up two particles moving on the branches of one definite 
hyperbola and at the same time have the particles fixed at the two available 
foci. 

There are three simple ways to view a binary system. Two of them 
are from the rest mass frame of one or other of the pair and the third is 
from the centre of gravity frame for the pair. This observation is valid 
for all combinations of gravitationally characterised components. I have 
chosen to carry through the calculation here within the rest mass frame 
of the positively gravitating mass m\. As earlier, I shall call vn\ the red 
component of the binary pair. In the rest frame of the red particle, the 
negatively gravitating green particle will be assumed to move on an 
hyperbola which as a whole is at rest, The path on which the green particle 
moves will be called the green path. This path will be an hyperbolic branch 
containing the stationary red particle fixed at the focus of that branch. The 
second fixed branch of this hyperbola does not play an active role so its 
focus will be displayed as an empty green circle in the diagrams 1 on page 
24 and 3 on pages 25 to be found at QMUL Maths, It will be helpful to 
refer to these diagrams for the following discussion. Working in the rest 
frame of the red particle, which is by definition fixed in this frame, will 
have the consequence that the path of the red particle, the red path, will 
have to be in motion, rotating and translating while changing points on the 
red path remain attached to the definite fixed position of the red particle. 
Thus for the frame we are working in, the red path moves through the fixed 
red particle position rather than the red particle moves on the red path. 
The green particle is on the green path but if it is to exert a repulsive force 
on the red particle it must also be at a remote focus of the red path. Thus 
mathematically we have to set up the two conditions: 

1 Red particle on local focus of green path 

2 Green particle on remote focus of red path 
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The hyperbolic parameterisations for various hyperbolic branches with 
respect to the two foci posibilities are given below, 





9,a,e) = a(e — cosh(#)) 


(8.19) 


xrl{ 


9,a,e) = — a(e + cosh(0)) 


(8.20) 


%lr{ 


9,a,e) = a(e + cosh(#)) 


(8.21) 


Xrr{9 


, a, e) = — a(e — cosh(#)) 


(8.22) 


Vll{9 


,a,e) = a{e 2 - l) 1/2 sinh(fl) 


(8.23) 




= VLR(9,a,e) 


(8.24) 




= VRL(9,a,e) 


(8.25) 




= VRR(9,a,e). 


(8.26) 



The two letter subscripts LL etc refer to the focus and branch involved 
respectively. For example LR means left focus for origin of coordinates and 
right branch for path. This terminology is OK provided we rethink when 
a rotating hyperbola turns through more than ±7r/2. Thus to impose the 
conditions (1) and (2) into the mathematics of our binary system, we note 
that the separation distance r s , between our red and our green particle can 
be expressed in two ways, in terms of the red particle parameters or in terms 
of the green particle parameters, each way giving the same numerical value 
for r s We are working in the red rest frame for the red particle and this is 
also the rest frame for the whole hyperbolic path of the green particle. This 
path will also have a fixed axis. However the red particle's hyperbolic path 
will be rotating in this frame so that we have to introduce its axial rotation 
away from the fixed axis of the green particle path. This rotation I shall call 
the angle (3. (3 is the angle between the green path's axis and the red path's 
axis. The positions of the various parameter components are given by the 
blue lines in diagram 3 on page 25 to be found at QMUL Maths together 
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with identification positions A, B, C, D, E, F denned below 



I? A 

t a 




xrl 


{O.Z I J 


FD 




Url 


(8.28) 


ZBFA 




P 


(8.29) 


ZEFD 


= 


P 


(8.30) 


FB 




x RL cos(/3) 


(8.31) 


ED 




Url sin(/3) 


(8.32) 


BA 




xr L sin(/3) 


(8.33) 


CB 




Vrl cos(/3) 


(8.34) 


DC 




x L l 


(8.35) 


CA 




Vll 


(8.36) 


FB - 


ED 


= DC 


(8.37) 


CB + 


BA 


= CA. 


(8.38) 



Thus using (I8.37j) and f!8.38j) and equating components for r s in the two 
parameterisations, we obtain 

XR L (0 R ,aR,eR)sm((3) + y RL {0 R , a R , e R ) cos(/3) 

= VLL(0G,a G ,e G ) (8.39) 

-Vrl(0r, a R , e R ) sin(/3) + x RL (9 R , a R , e R ) cos(/3) 

= x LL (9 G ,a G ,e G ). (8.40) 

My objective is to produce a single nontrivial case of a mixed gravity bi- 
nary system to establish that such systems can theoretically exist without 
internal contradiction or violation of physical laws. This can be most easily 
achieved by taking the simplest case. We have seen that the masses of the 
two components need be different and this translates into the a's being not 
equal. The other parameter e can be taken equal for the two orbits just to 
reduce notation and the arithmetic. Thus from now on I shall take a G = a 
and a R = 0.3a G = 0.3a. Then the basic equations become 

xrl(0r, 0.3a, e) sin(/3) + Vrl{6 r , 0.3a, e) cos(/3) 

= y L L(d G ,a,e) (8.41) 

-Vrl{Qr, 0.3a, e) sin(/3) + x RL (9 R , 0.3a, e) cos(/3) 

= x LL (9 G ,a,e). (8.42) 
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These last two simultaneous equations can be solved for the pair of variables 
sin(/3), cos(/3) to give 

sin /m xrl(0r, 0-3a, e)y LL (9 G , a, e) - y RL (9 R , 0-3a, e)x LL (9 G , a, e) 

xrl{0r, 0.3a, e) 2 + y RL {9 R , 0.3a, e) 2 

(8.43) 

Vrl(0r, 0.3a, e)y LL (6 G , a, e) + x RL (9 R , 0.3a, e)x LL (9 G , a, e) 



cos(/3) 



%rl{6r, 0.3a, e) 2 + y RL (9 R , 0.3a, e) 2 

(8.44) 



The expanded versions for sin(/3) = S(9 R ,9 G ) and cos(/3) = C(9 R ,9 G ) are 
given next in terms of the parameters, 9 R and for position on the two 
hyperbolae, 



sin OS) 



e 2 - iy/ 2 \ /cosh(fl G ) smh(0 R ) - sinh(fl G ) cosh(^) 



0.3 J V (ecosh(fl H ) + I) 2 

e(sinh(6» G ) + sinh(^)) 
(ecosh(^) + l) 2 
S(9 R ,9 G ). (8.45) 



cQ 3 /e 2 (sinh(^) sinh(fl G ) - 1) - e(cosh(0 fl ) - cosh(£ G )) 



10 V. (ecosh(^) + l) 2 

sinh(#R) smh(9 G ) + cosh(6'^) cos1i(#g) 



(ecosh(0 fl ) + l) 2 
= C(9 R ,9 G ). (8.46) 

Thus it is possible to find the relation between 9 R and 9 G for the binary 
system to hold together consistently by using the identity, 

sin 2 (/3) + cos 2 (/3) = 1 (8.47) 

and employing the functions above at (I8.45P and (18.461) . However, the for- 
mula that emerges by that route is intrinsic and complicated in the variables 
9 R and 9 G and it is difficult to see how to extract an explicit relation giving 
for example a function of one in terms of the other such as, 9 R (9 G ). A useful 
result that can easily be obtained from the S and C functions is the angle 
between the two orbit axes for any pair of positions. This is given by the 
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inverse tangent 



^'"^-'(ra) (8 - 48) 

and can be used in finding the moving rotating path of the red particle in the 
green frame or vice versa. If we return to the original pair of equations, we 
can find the relation between the two position on curve parameters, 9 R and 
9 G , simply by squaring the two equations and adding the results together 
to obtain, 

rloi^R, e ) = x RL (9 R ,a R ,e R ) 2 + y RL (9 R ,a R ,e R ) 2 (8.49) 
= x LL (9 G ,a G ,e G ) 2 + y LL (9 G ,a G ,e G ) 2 (8.50) 
= rl R (9 G ,e), (8.51) 

whilst eliminating /3 at the same time. This result simultaneously defines 
the square of the separation distance, r 2 G (9 R ), between the two particles, 
when viewed from the green particle or as, r 2 R (9 G ), when viewed from 
the red particle, 9 R varying from the green point of view and vice versa. 
Expanding the expressions in (18.491) and (18.501) we find, 



rl G {0R 


e) 


= (0.3a) 2 (ecosh(#R) + l) 2 


(8.52) 


r 2 s , R (e G 


e) 


= a 2 (ecosh(# G ) - l) 2 


(8.53) 




e) 


= ±0.3a(ecosh(fl R ) + 1) 


(8.54) 




e) 


= ±a(e cosh^c) — 1) 


(8.55) 




>g) 


= ± cosh~ x (± cosh(fl G )/0.3 T l/(0.3e) - 1/e) 


(8.56) 




»h) 


= ± cosh _1 (±0.3 cosh(^) + 1/e ± 0.3/e). 


(8.57) 



The fifth equation above is the sought relation between 9 R and 9 G resulting 
from the equality of the two separation views from the first two equations. 
The leading ± at (18. 561) and (18.571) results from the evenness of the cosh 
function with the ~ indicating that this ± is independent of the others, 
so implying four possibilities. The sixth expression above is the inverse 
transformation of the fifth. The polar forms for the separation distances 
can be obtained by replacing the parameters 9 R and 9 G by functions of 
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their angular equivalents 4> R an d 4>g, say. 



a (a \ *. -1 f-ynL(0 R ,0.3a,e)\ 

(j) R {0 R ) = tan 8.58 

V x iJL (6' R ,0.3a,e) / 

a Cfl ^ + -i f VRL{0 G ,a,e)\ 

= tan U M (fe.«.«) J ( ' 

= to,- ( ""<*■'■'■'» )■ (8.60) 
\XLL{QG,a,e) J 

Both expressions at ( 18.5411 and ( 18. 55ft after the ± signs above are always 
positive because e > 1 and cosh(a;) > 1. Consequently, one might infer that 
the negative options can be ignored on the grounds of being unphysical, if 
the r SjR /Q quantities represent physical distance. This turns out to be not 
the correct inference so that the final definitions of these quantities have 
to be such that they, being scalar distance, are carefully defined to come 
out positive. This will be explained. The inversions of the three equations 
above, the 6s in terms of the 0s are 

a (a \ w /etan 2 (0 fi ) - sec(0ij)(e 2 - 1)\ 

M« = cosh ^ e , _ secH4>R) ) < 8 - 61 > 

fete) = cosh-'f ^^'^+^f- 1 ' ) (8.62) 
V e 2 -sec 2 (0 G ) / 

„. , i /etan 2 (0) — sec(0)(e 2 — 1)\ . 

W = «"*-'( el-JZ) "J'' <8 ' 63) 

As a result of these equations, we can re-express the formulae for the dis- 
tance between red and green particles, the r StR (9c) and r St c(0 R ) functions 
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at f)8.54p and (I8.55p . in terms of the angles 0g and 4>r respectively. 



— >■ 



±O.3a(ecosh(0 fl (<fo)) + 1) 

e 2 tan 2 (0 jR ) — sec(0 jR )e(e 2 — 1) 



0.3a 



e 2 - sec 2 (0 jR ) 



0.3a(e 2 - r 



r U(0G,e) 



cos(^ R )e + 1 
±a(e cosh(0 G ) — 1) 

e 2 tan 2 (0 G ) + sec(</>G)e(e 2 — 1] 



e 2 — sec 2 (0c) 



a(e 2 - 1) 



cos(0c) e — 1 



8.64) 
8.65) 

8.66) 
8.67) 
8.68) 

8.69) 



the primes being introduced to differentiate these new functions of (f) from 
the original functions of 9. The ± options have been assigned so that both 
these length functions are positive as indicated by the modulus sign, | 
enclosing their final definitions. This is the point I suggested earlier should 
be treated with care as the introduction of the versions with the inverse 
coshes seems to open up the possibility for negative values. The relations 
between the 4>rs and 0gs corresponding to the relations between the 9rs 
and 9 G s at (I8.56P and (I8.57P which make the two lengths equal are obtained 
easily by taking the two lengths as equal and then solving for either <pG or 
4>r and are 



Vg( 
M 



= cos 



1 ((cos(fo)/0.3 + 1.3/(0.3e))) 
= cos _1 ((O.3cos(0 G ) -1.3/e)) 



(8.70) 
(8.71) 



This completes most of the technicalities and leaves us with the two dis- 
tances, r' s g (0_r) and r' s r (4>g), between the red and green particles as would 
be seen by an observer on the green particle in terms of the red particle's 
angular parameter <J)r and an observer on the red particle in terms of the 
green particle's angular parameter (fie respectively. The green observer sees 
the red particle and its path and the red observer sees the green particle and 
its path. However, these distances are not of equal magnitude unless the 
basic Or and the basic Qq are related by the formula (I8.56P or the formula 
(18.571) or the equivalent 4>r and (fie are equivalently related. 
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9 Assembling the Paths 



Most of the mathematics for this problem has now been completed. We 
have seen that in the rest frame of the red particle, the path of the green 
particle is also at rest, the red particle being at the focus of this hyperbolic 
path. It is also established that the axis of the path of the red particle will 
have to have rotated relative to the green path axis by the angle (5 given at 
(I8.48p . Concentrating on the view of the situation from the rest frame of the 
red particle in which the path of the green particle is a simple hyperbolic 
branch also at rest, it is convenient to introduce a rotated frame, relative 
to this at the angle /5 to it, in which the green particle is at rest and the 
path of the red particle is also at rest. To do this, I introduce the rotation- 
translation transformation of coordinates based at the empty green focus 
which is at (0,0), 



where d\ and d\ are displacements following the rotation that take the empty 
green focus to the momentary position of the green particle assumed to be at 
parameter value Qq- Thus this transformation takes us into a reference frame 
in which the green particle is at the origin of the transformed x-axis which 
is at the angle /3 relative to the axis of the green parabolic path. In other 
words, the x-axis of this frame is on the x-axis of the red hyperbola with the 
green particle on its active focus. Thus to find the parametric equation for 
the red particle path all we have to do is place its parametric representation 
into the rotations (x, y) coordinate position and replace the rot subscripts 
with red to indicate we now have the two parametric components for the 
path of one of the branches of the moving red hyperbola as functions of (f)R 



X rot (x,y,d 1 ,f3) 
Y rot (x,y,d 2 ,(3) 
d\ 
d% 



(xcos(/3) — ysin(/3)) + d\ 
(xsin(/3) + ycos(/3)) + d 2 

y RL (0 G ,a,e), 



(9.1) 
(9.2) 
(9.3) 
(9.4) 
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and 9 G . 



X re d{<p R ,9 G ) = x RL (9 R (4> R ),0.3a,e)cos(fi(9 R ((j) R ),9 G )) 
y RL (0 R ^ R ),0.3a,e)sm((3(9 R ( ( p R ),9 G ))) 

+ x RL (9 G ,a,e) (9.5) 

Y red ((f> R ,9 G ) = x RL (9 R ((f) R ),Q.3a,e)sm([3( y 9 R ((j) R ),9 G )) 
+ y RL {9 R (4> R ), 0.3a, e) cos(P(9 R (4> R ), 9 G ))) 

+ yRL(9 G ,a,e), (9.6) 

where 9 R has been related to its angular equivalent (fi R and (3 has been 
replaced by the function of 9 R and 9 G at (I8.48P repeated below. 

The other co-moving branch is given by 

Kedi^R^a) = x RR (9 R ((f) R ),0.3a,e)cos(l3(9 R ((f) R ),9 G )) 

- y R R{9 R {<f)R), 0.3a, e) sm(P(9 R (<f) R ), 9 G ))) 

+ x RL (9 G ,a,e) (9.8) 

KMr^g) = x RR (9 R ((f> R ),0.3a,e)sm^(9 R ( ( f> R ),9 G )) 

+ yRR(0 R {4> R ), 0.3a, e) cos((3(9 R (<p R ), 9 G ))) 

+ y R L{9 G ,a,e), (9.9) 

with the same adaptations as for the first branch. With the parameterisation 
of the moving hyperbolic path for the red particle obtained at H9 .51) and 
(19. 6p . it is possible to plot that path which has the moving green particle 
at its active focus and so give a clear diagrammatic picture for any chosen 
position, the value of 9 G , of the green particle on it own path which you will 
recall is fixed, together with the red particle at its left focus, in the basic 
reference frame on which we are concentrating. The fixed path branches for 
the green particle is given by the parameterisations, 

X green 

(<j) G ) = x RL (9 G ((f) G ),a,e) (9.10) 

Y green ((j) G ) = y R L{.9 G (<p G ) } a } e) (9.11) 

KeMa) = x RR (9 G (<p G ),a,e) (9.12) 

Y'green^G) = y RR (0 G (<j> G ) , a, e) . (9.13) 
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We now have all the basic mathematical structure to plot and analyse the 
path structure and the relative motions of the particles in terms of the an- 
gular changes of <pR and or <f>Q. Most of the emphasis here has been the 
view of the binary pair of an observer fixed to the rest frame of the red 
particle and the accompanying fixed hyperbolic path of the green particle. 
From that frame of reference the observes sees only the motion of the green 
particle around its orbit. It is not difficult to recast the whole structure in 
terms of the view of an observer on the green particle and his accompanying 
fixed hyperbolic orbit of the red particle. Further these two basic views can 
then be put together using the definitions of the centre of mass or centre of 
gravity systems to give an observers view from either the centre of mass or 
the centre of gravity frame. However, only the last of these reference frame 
possibilities could be taken to be truly inertial. Using the Mathematica ani- 
mation process, I have produced a simulation of the movement of the mixed 
character binary pair with 301 frames. One of these frames has been used to 
produce the diagram on page 24. The full simulation in mathematica note 
book language can be downloaded in the file mixmass.nb from my website 



at, QMUL Maths| It can now be shown that the structure of the force 



configurations for this binary pair composed of oppositely gravitationally 
characterised particles does conform to the correct conditions. They are 
that the red particle exerts a gravitational attraction on the green particle 
and the green particle exerts a gravitational repulsion on the red particle. 
This can be done by using a well know result from classical dynamics and 
inverse square law forces which recovers the force involved from just know- 
ing the distance in detail of a possibly influenced particle from the force 
source as a function of the angle (p. The formula required is derivable from 
(HD and is 

a r>R = r ■ d 2 r/dt 2 = (f - rcj) 2 ) = -Gm% /r 2 (t) 

= h 2 u 2 {<f)){d 2 u{<j))/d<j) 2 + u{<f))) = -Gm+u 2 {<f)) (9.14) 



u(4>) = l/r(t) (9.15) 
h = r 2 (t)0 2 , (9.16) 

where normally m^, as indicated, the gravitating source would be posi- 
tively gravitating, h is the constant arising from integrating the transverse 
component of acceleration on the assumption that the mutual interac- 
tion of the particle pair only acts along their line of separation. As we are 
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working with the inverse square law of gravitation u 2 (<f)) appears twice in 
the formula (I9.14p and so can be cancelled from both sides of the second 
equality to define G t y P e(4>) and to give, 

G tyP e(<j)) = h 2 (d 2 u((j))/d<j) 2 + u(<f>)) = -Gm+ } (9.17) 

where u(4>) is the inverse of the magnitude of a position vector for distant 
points originating at the position of the mass m\. This shows immediately, 
because of the —G, that we are dealing with an attraction towards the 
source under gravity which in the past has been thought to always hold. 
The detailed steps in deriving the formula (I9.17P outlined above can be 
found in A. S. Ramsey's book Q42J). They are elementary but I think this 
formula needs to be treated with care because it is rather profound. In 
particular, it should be noted that the force involved on the right hand side 
above is the force originating from the particle 777,2 which is at the origin of 
coordinate from which the radial length, r(<p) = 1/m (</>), is measured and 
which is also the magnitude of a vector, r(0), which has its tail at the same 
origin. However, the angular parameter <p on which it depends refers to the 
angular component of acceleration at the sharp end of the same vector and 
is involved with the path where the induced distant acceleration field may 
be influencing other particles but even so this (f) is an angular parameter 
measuring an angle as perceived by an observer at the source of the force. 
According to general relativity all such other particles at the sharp end of 
the vector r(<p) will experience the same acceleration due to the gravitational 
field of the particle 777.2. So there is some sense that general relativity is at 
least a part feature of the formula. Thus when using this formula to find the 
nature of the force involved we must not lose sight that it refers exclusively 
to the source location of the force. The formula determines how the source 
observer sees the process at a distance but it contains no feedback from 
the distant events. This has the consequence that the function u{<p) which 
is used in the formula should be specific to the source as origin which one 
is trying to interpret. That is to say, it should not be in a form which 
makes it equal to the distance between the binary pair as observed from 
the other member of the pair. If the physically equal distances are made 
mathematically equal before using the formula, The two u(<f>) functions 
become identically equal and so the differential form in the expression (19.171) 
gives the same result for the two binary members. The significance of this 
formula is that it gives information about the local source force, whether 
attractive or repulsive, valid for all receptor particles at whatever distance 



34 



from the source they may be and whatever gravitational character they may 
have. This independence on receptor distance is indicated in the diagrams 1 
and 2 on page 25 to be found at QMUL Maths by plotting the formula value 
with corresponding value of Thus we can use the raw functions of 4> R 
and (fic at (18.66j) and ( I8.69P to determine the type of gravity originating 
from the green particle or the red particle respectively These expressions 
are repeated below 



r a ,G\.<PR> e ) 



0.3a(e 2 - 1} 



cosO/j)e + 1 
a(e 2 - 1) 



cos(0c) e — 1 



l/u Green ((f) R ,e) (9.18) 
l/u Red ((P G ,e), (9.19) 



with the second equalities giving the appropriate u((f), e) function. If we 
evaluate the first and second derivatives of UGreen((f)R,e) and UR e d(<j>G, e), 
with respect to 4> R and 4>g respectively and then use them to find the value 
given by G type for the two cases, we find 



Gfype ( 
Gtype ( 



'>R) 



h 2 (d 2 u Gree 
h 2 (d 2 u Red { 



^R, e )/d(j) R + ^Green 



e)/d(p 2 G + u Red ( 



R,e))>Q 
G,e)) <0, 



(9.20) 
(9.21) 



implying that Gt ype {(f>R) = Gm + and Gt ype {(f>G) = —Gm~ and thus show- 
ing, as expected, that the red particle will attract all other particles in its 
neighbourhood and the green particle will repel all other particles in its 
neighbourhood. The situation can be clearly seen from plots of G t y P e((pR) 
for the two cases that are given at diagrams 1 and 2 on page 25 to be found 
at QMUL Maths. The corresponding u((f>) for the two cases is plotted with 
the G tyP e{(t>) block function and this shows clearly how the character of any 
source does not depend on r(<f>) = l/u(<p), the distance of any other particle 
in its vicinity whatever character that particle may have. 



10 Conclusions 

The work in this paper makes strong theoretical evidence for the possible 
existence of particulate negative gravitating positive mass particles. This 
has been achieved by showing that negatively gravitating particles can be 
involved in composite systems and follow orbits in conjunction with posi- 
tive gravitating particles in much the same way as the later type of usual 
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particles can form systems and gravitationally interact. In the specific case 
examined here, this new orbital system type does not have the bizarre fea- 
tures that were noted in the fifties studies by Herman Bondi. In the light 
of this special but typical case, it seems reasonable to assume that if such 
particles do exist their interaction generally with the normal type of parti- 
cles would be theoretically predictable and not remarkably odd. Of course, 
I have not given any physical evidence that such particle do exist. 

There are a number of interesting issues that arise from this work. One 
such issue is the meaning and significance of the term particle. We talk 
about particle in quantum theory and at the other end of the scale, in 
cosmology, galaxies are often regarded as particles. In the definition of 
particle that has been used in this work the whole of this vast range has been 
included because it is generally believed that all such particles whether they 
have rest mass or not are influenced by and can cause gravity fields. However 
we think we know that negative gravity particles, if they exist at all, will be 
very rare entities. It is a fact that antiparticles in general are very rare, a fact 
that is regarded as a great mystery. The question is, where are the missing 
antiparticle? Some scientists believe that normal and anti particle were 
created together in equal numbers at an early stage of the evolution of the 
universe but now most of the antiparticle seem to have disappeared. This 
prompts me to make the speculation that dark energy is composed of those 
missing antiparticles and further those antiparticles that we do detect are 
manifestations of the local sea of dark energy particles and like them are also 
anti gravitational. This view is also supported by the accepted cosmological 
fact that at some time t c in the evolution of the universe there was equal 
quantities of normal mass and dark energy mass within the boundaries of 
the universe. This time could be the moment of creation of the equal anti 
and normal particle densities. The same identificaton for dark energy was 
put forward as a conjecture by D. S. Hajdukovic in October 2008 [33]. 

Staying with the particle concept problem, there is second interesting 
issue. How is it that accumulations of mass arise? It is a convincing idea in 
cosmology that large accumulations of mass such as is found in a in a galaxy 
could have arisen by gravitational accretion or by a dispersed distribution 
of particulate mass falling together as a result of the mutual gravitational 
attraction of its component particles. Stars could be formed in the same 
way. This may be what happens in some or many cases but it is possible 
that such structure are formed and spill out from some unexplained minute 
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but greatly massive singularity. Whatever the true situation is, it would 
seem that a dispersed distribution of negative gravitating particles would 
not accrete into a single mass because the elements of the distribution would 
be mutually repulsive. This suggests that negative gravity galaxies of the 
usual type are unlikely to exist. However, there is another possible twist 
to this issue and that is how do the mass accumulations within elementary 
particles form. Some researchers think that gravity is involved in the struc- 
ture of elementary particle. They could have been formed by some sort of 
gravitational collapse and it could be that gravity holds them together. The 
snag is that gravity is a very weak force and such ideas are not at all con- 
vincing. In fact, we have no idea what holds elementary particle together or 
how they are first formed from energy. In quantum, mechanics however, we 
do have good ideas about how particle transform between themselves. This 
does have a bearing on the question of possible galactic size negative mass 
particles. Dense matter that is found in atoms, molecules, rocks, planets, 
stars and galaxies is held together by atomic forces of one sort or another. 
The smaller accumulated parts then held together in the large by positive 
gravity. This suggests that large accumulations of negative mass particles 
could be possible if constructed by atomic processes without the interven- 
ing spaces found in normal galaxies. I am suggesting that galactic sized 
molecules could be formed from negatively gravitating particles. Enough 
speculation, I return to the main issue of the orbiting binary pair and its 
lessons on fundamentals in the next paragraph. 

Another issue of interest is the question of Newton's Action equals Reac- 
tion principle, his third Law, and how that squares with Action at a distance 
in the case of the gravitational interaction between separated positive and 
negatively characterised particles. Other than mentioning this issue, I have 
carried through the work in this paper as though there is no problem. In 
fact, that is the case, there is no problem this being an important con- 
sequence that follows through from Einstein's principle of equivalence in 
general relativity. 
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To explain this issue, I first quote one popular version of Newton's third 
law: 

// a first body exerts a force f on a second body, 
the second body exerts a force -f on the first body. 
f and -f are equal in size and opposite in direction. 

Alternatively, in a form nearer to Newton's original presentation: 

All forces occur in pairs, and these two forces are 
equal in magnitude and opposite in direction. 

Newton's original version of this law is not very specific in respect of what a 
force is and what it is doing but it does seem to be talking about one basic 
force and its somehow reflected mate such as we all experience when we push 
an object with our hand and feel the resisting response, presumably the 
reaction to our action. Such a situation is decidedly local and is essentially 
all about a single basic force and is very intuitively convincing. The first 
definition talks about bodies, says essentially the same as does Newton's law 
but although it mentions bodies it does not define a body or explain how 
these two bodies are orientated in space. Of course, fundamental definitions 
have to be minimal. Neither of these definition obviously apply to bodies 
separated in space in the way that a binary pair of particle, whatever their 
gravitational characteristics may be, are separated in space. However, it is 
usual practice to invoke the action and reaction principle in the standard 
binary pair case to justify saying the force on one from the other must be of 
opposite sign to the force from the other to the one. This practice does give 
a correct result in the usual binary pair case but as we have seen it seems 
not to work for the mixed mass pair. This situation prompts me to suggest 
that calling on Newton's third law to support the proposition that the two 
possible gravitational forces between distant particle from one or the other 
must be of opposite sign is philosophically incorrect, in spite of it giving a 
correct result. Before I make this case, let us consider how general relativity 
theory describes the properties of and incorporates a local gravitational 
field arising from a distant source. Local acceleration fields generated by 
distant gravitating particles take precedence over the gravitational forces 
experienced by any local particle. This is essentially Einstein's Principle of 
equivalence which manifests itself by the fact that all particles in a given 
location move with the same acceleration due to local gravity, regardless of 
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their individual masses and, indeed, this is also part of the Newtonian theory 
structure as described by equations (17.81) and (17.121) from which it is clear 
that the mass of a gravitationally influenced particle cancels out from both 
side of the equation of motion under gravity. Essentially this means that two 
spatially separated gravitating particles produce independent acceleration 
fields in which the other particle moves. Each of the separated particles 
causes an action at a distance on the other particle but the two particles 
are acting independently. This is certainly not reaction of one to the action 
of the other as implied by Newton's third law from whichever of the two 
possible directions it may be viewed. However, I have shown above that in 
the classical case the forces involved do satisfy the condition 

fi = -f 2 , (10.1) 

which is the condition for action and reaction to be equal in this classical 
case of a positively characterised pair. Thus if the equal action reaction 
property holds, how is it reconcilable with my claim that this is not an 
example of Newton's third law? The answer to this question is that there is 
a more fundamental simple and kinematic reason of greater generality than 
Newton's third law. We can see this clearly by examining the equations at 
(I7.32p and (I7.34p rewritten below in a rearranged form. 



+ d 2 r(t)/dt 2 = 


fi f 2 


(10.2) 


M+ M+ 


M + = 


mfmf 


(10.3) 


mf + mf 


-d 2 r(t)/dt 2 = 


fi f 2 


(10.4) 


M_ M 


M_ = 


mfmf 


(10.5) 


mf — mj 



where M + and M_ are the appropriate reduced centre of mass or centre of 
gravity masses for the system under consideration. 
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The relative separation acceleration vectors for a binary system as seen 
from each particle of the pair for the classical case, the first two, and for 
the mixed mass case, the second two, can be defined as, 

a 1+ = A- (10.6) 

a 2+ = A (10.7) 

a,. = A- (10.8) 

a 2 - = A (10 .9) 

If negatively gravitating particles do exist then, according to (110.21) and 
(I10.4p all the situations that can occur with a same type gravitating pair or 
a mixed pair can be summarised as, 

(10.10) 

The formula (110. 10p for relative accelerations has a very simple kinematic 
and indeed obvious explanation in terms of the geometrical distance between 
the particle pair, their centre of gravity or their centre of mass. In the case 
of a same characterised pair, the unit vectors ri and f 2 directed from the 
members of the pair towards their centre of mass which is between the end 
points of pairs separation distance and so must be of opposite sign ri = — r 2 . 
This leads to the fi = — f 2 and then implies, with the usual terminology, 
Newton's action is minus reaction condition. However, in the first case 
the terminology action equals minus reaction is a misnomer so that the 
second case when fi = f 2 is not in conflict with Newton's third law. In 
the second case of a differently characterised pair, the unit vectors ±f 1 and 
±r 2 directed from the members of the pair towards their centre of gravity 
which is outside the pairs separation distance and so must be of the same 
sign ±ri = ±f 2 . This leads to the fi = f 2 and then only seems to imply the 
violation of Newton's actions is minus reaction condition of the first case. 
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I suggest that in astronomical physics in the context of separated bodies 
Newton's law of action and reaction could usefully and more realistically be 
replaced with the new law:- 

The relative gravitational separation acceleration vectors ai and &2 
for two bodies satisfies one or other of the two case, ±, relation 

a 1± = T a 2 ±. (10.11) 

Whether negatively gravitating particles exist or not Newton's law of ac- 
tion and reaction in the astrophysics of separated bodies could usefully be 
replaced as above with only the + sign in the subcript and the minus sign 
multiplying. 

I have carried through the work in this paper using only Newton's equa- 
tions of motion but with much guiding influence from general relativity 
ideas. I think a fully relativistic treatment of this topic would be very com- 
plicated and not reveal any substantial deviations from the conclusions that 
I have arrived at here. 
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